A new perspective on the Holstein polaron problem 
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The single-polaron band structure of the Holstein model in one and two dimensions is studied 
using a new form of resummed strong-coupling perturbation theory. Well converged results are 
obtained for phonon frequencies of the order of the hopping integral and strong to intermediate 
electron-phonon coupling. The polaron band structure at intermediate coupling is shown to deviate 
markedly from that of a nearest-neighbor tight-binding model, and is in fact similar in shape to the 
prediction of weak-coupling self-consistent perturbation theory. 
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Despite many years of research, our understanding of 
the nature of the transition from a quasi-free electron 
state with small mass renormahzation at weak electron- 
phonon (EP) coupling to the small polaronic state with 
very narrow band width at strong coupling is still incom- 
plete. In particular, the region of intermediate coupling 
and phonon frequency is difficult to handle theoretically 
due to the absence of a small parameter on which to base 
a perturbation theory. It is known rigorously that the 
ground state energy and effective mass in the Holstein 
modetl are analytic functions of the EP couplingcl, which 
implies that there is no true phase transition in this sys- 
tem. Nevertheless, accurate Monte Carlo calculations^! as 
well as direct diagonalization approaches^ demonstrate 
that the character of the ground state changes signifi- 
cantly at a rather well-defined crossover coupling. 

In this paper we concern ourselves with the evolution of 
the polaron quasiparticle (QP) band structure as a func- 
tion of EP coupling. We present the results of a study of 
the Holstein model using a new numerical scheme, which 
is equivalent to a selected resummation of some terms 
in strong-coupling perturbation theory (SCPT) to all or- 
ders. Our resummed SCPT, denoted FC for "finite clus- 
ter" , since the method is a form of cluster expansion, con- 
verges well down to quite weak coupling also for phonon 
frequencies of the order of the electronic hopping inte- 
gral, where we are able to obtain good agreement with 
the weak-coupling self-consistent Migdal (SCM) pertur- 
bation theory and exact results for small clusters of up to 
six sites. We present a simple physical picture which clar- 
ifies the continuous evolution of the QP band structure 
from weak to strong coupling. 

The Holstein modeltl is defined by 



H = -t ^ (c}Cj + C^Ci) +UJ ^ a i a i 
<i,j> i 

+ gJ2m(4+ a i) (!) 



where cj is the creation operator for a spinless electron 

on site i,m = cfc, and a\ is the creation operator for the 
local oscillator of frequency u>q. Although the formalism 



presented may also be applied to the many-electron case, 
we treat only the single carrier case here. 

We recall first some well known results for the strong- 
coupling limit: if the EP coupling is large it is natural to 
treat the coupling term in (|l|) exactly, and to consider the. 
hopping as a perturbation. Performing the Lang-FirsovO 
canonical transformation leads to H' = Hq + H t , where 



(2) 



is diagonalized by the transformation, and the kinetic 
energy term becomes 



H t = -t 



<i,3> 



(3) 



with Xi = exp a{a,i — <xj) and a — g/ojQ. 

One now chooses as the model space states with no ex- 
cited phonons, and upon taking the expectation value of 
(^) in the transformed phonon vacuum (first order SCPT) 
arrives at an effective polaron Hamiltonian consisting of 
a nearest-neighbor (n.n.) tight-binding model with the 
hopping integral t — > t* = <cxp(— a 2 ) together with a 
constant energy shift (binding energy) of — Eb — —a 2 u>o. 
Continuing withm-degenerate Rayleigh-Schrodinger per- 
turbation theoryQ'EL the second order correction involves 
the application of (^|) onto an initial state with no excited 
phonons, leading to an intermediate state with the elec- 
tron on a n.n. of the initial site and some excited phonons 
as well. The second application of (||) must annihilate the 
phonons and either a) return the electron to the initial 
site, leading to a renormahzation of the binding energy or 
b), move the electron to a next-nearest-neighbor (n.n.n.) 
of the initial site, giving rise to an effective n.n.n. hop- 
ping. Making use of translational symmetry one finds in 
one-dimension (ID) the band structure 



E(k) = -E B + E {2) - 2t*cos(k) - 2t 2 cos(2k) (4) 



where 



E^ = 



-2{t* 

UJQ 



[Ei(2a 2 ) - 7 - ln(2a 2 )] 



(5) 



1 



«2 = — [Ei(a 2 ) - 7 - ln(« 2 )] , (6) 

with Ei(x) the exponential integral and 7 Euler's con- 
stant. 

Note once more that the initial tight-binding model 
contains only n.n. hopping terms, but due to the strong 
coupling to the lattice degrees of freedom the polaron 
dispersion develops n.n.n. termsEl . At higher order in 
perturbation theory one generates further longer range 
hopping terms as well as renormalizations of the previ- 
ously obtained contributions. If one considers uio s=s t the 
convergence of this series is poor, even for quite strong 
coupling. One option for a systematic and controlled 
approach to a description of this intermediate regime is 
to continue as we have started and to generate further 
terms in the Rayleigh-Schrodinger perturbation expan- 
sion. While this is clearly feasible, it may be necessary 
to go to quite high order to achieve converged results for 
a wide range of parameters. We choose instead an alter- 
native approach: it has been shownO that it is possible 
to resum to infinite order some of the terms of a strong- 
coupling perturbation expansion by making appropriate 
use of some exact information obtained from the numer- 
ical diagonalization of a sequence of finite clusters. This 
method has been shown to allow the low-lying states of 
the Hubbard model to be quantitatively well described 
by a generalized t — J model also for Uf t not very largem. 

The general form of the result obtained is a series of 
terms involving connected clusters of increasing numbers 
of sites, h = + + + where is an op- 
erator which acts on all connected clusters of exactly n 
sites. By construction, the eigenvalues of h will be exact 
eigenvalues of the original full Hamiltonian on this same 
cluster if terms up to ftW are applied to a cluster of n 
sites. 

The terms of h are constructed recursively: given the 

series up to h^ n ' one may extract a "total" h for an n + 1 

site cluster of a particular topology as h = UEU -1 , 

where E is a matrix containing a subset of the exact 

eigenvalues of the full original Hamiltonian for this given 

cluster on the diagonal, and the transformation matrix 
1/2 

is U = (5~ 1 t<S'~ 1 ) S. Here S is the matrix whose 
columns are the projections of the relevant eigenvectors 
of the full Hamiltonian onto the model space (in the 
present case the zero phonon subspace). The factor of 

ensures that h is Hermitian (for fur- 
ther details and explanation see Ref. |9|). 

The exact information for E and S is obtained in the 
present case by numerical diagonalization after first ren- 
dering the problem-finite by introducing a cutoff in the 
phonon occupation^ The matrix elements of are 
now found from those of this "total" h on the (n+ l)-site 
cluster by subtracting the contributions from all embed- 
ded subgraphs of up to n sites which have previously been 
calculated. It should be stressed that this procedure is 
not a fitting or extrapolation of eigenvalues: in general 



also information from the exact eigenvectop,enters in the 
determination of the matrix elements of MlH 

For the single carrier case the resulting effective Hamil- 
tonian is a single-particle tight-binding model which may 
be diagonalized directly for the infinite system. If the 
same approach is applied to a many-carrier system one 
may "integrate out" the phonon degrees of freedom just 
as in the single particle case, but one is then still left with 
a complicated many-fermion problem to analyze. Work 
in this direction using numerical diagonalization of the 
resulting effective Hamiltonians for finite systems is in 
progress. 

The convergence of the SCPT for the Holstein model 
which is evaluated up to second order in (Q) depends on 
two parameters: both the degree of adiabaticity ui^/t and 
the coupling strength a are important. For tuo/t >> 1 
(anti-adiabatic limit) the expansion is well convergent for 
all a, but for u>o/t ~ 1 large a is required. This behav- 
ior is also reflected in the resummed FC approximation 
described above, where it manifests itself in the spatial 
range of interaction needed to achieve a converged re- 
sult. In general the smaller the phonon frequency the 
larger the interaction range needed for a given a. 

We compare in Fig. 1 the predictions of various ap- 
proximations to the polaron dispersion in ID for an inter- 
mediate phonon frequency of u>o/t = 1 and two different 
coupling strengths. 
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Fig. 1. Polaron dispersion for a ID system with ojo/t = 1 
and: (a) g/t = 1.75; (b) g/t = 0.8. Solid lines are FC ap- 
proximation to order /i' 5 ', dotted lines FC to h^K Dashed 
lines are second order strong-coupling perturbation theory, 
and dot-dashed lines weak-coupling self-consistent Migdal ap- 
proximation. Data points are exact results for various small 
clusters with periodic boundary conditions: circles, 6 sites; 
squares, 5 sites; triangles, 4 sites. 

The extremely small difference between the predic- 
tions of the 4- and 5-site FC approximations in Fig. 
1(a) demonstrates the rapid convergence of this approach 
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in strong-coupling even for this "difficult" phonon fre- 
quency. These curves are also in very good agreement 
with the diagonalization results, with the discrepancies 
being of the same order as the finite-size scatter from one 
system to the next. The second order strong-coupling 
perturbation theory on the other hand underestimates 
the band width by approximately a factor of two even 
for this relatively strong EP coupling because of the low 
phonon frequency. The weak-coupling SCM result, which 
is the lowest pole of the single-particle Green's function 

G(k,w + iS) = [uj + iS-e k -S(fc,w)] _1 (7) 

with the self-energy 

S(fc, W ) = SH= S 2 /^EG(fc>- W o) (8) 

k' 

is as expected quantitatively very unreliable here as well. 
MarsiglicQ has previously compared the ground state en- 
ergy predicted by this approximation with exact numeri- 
cal results, and found good agreement only for quite weak 
coupling. We have furthermore confirmed that adding 
the first vertex correction to the SCM0 does not signif- 
icantly improve the agreement of the band width with 
the results of the other methods for parameters as in 
Fig. 1(a). 

Turning now to the weaker EP coupling case shown in 
Fig. 1(b), we see that the FC predictions now agree quite 
well with the weak-coupling SCM prediction for the QP 
dispersion, as well as with the small-system data points. 
In this case the FC approach displays slightly poorer con- 
vergence than in Fig. 1(a), with the curves for the two 
different interaction ranges shown oscillating noticeably 
about one another. For these parameters the polaron 
size is slightly larger than the range of interaction in- 
cluded in the present calculations, but nevertheless it is 
clear that the "limiting" result would not be very dif- 
ferent from those shown. One expects predominantly a 
further smoothing of the oscillations with increasing ef- 
fective interaction range. Note once more that the FC 
approach starts from the strong-coupling limit, so that 
the good agreement with the dispersion predicted by the 
SCM weak-coupling approximation and the exact finite- 
system data in Fig. 1(b) is non-trivial. In this case the 
simple second order SCPT prediction is quite poor, but 
errs in the opposite direction from the stronger coupling 
case of Fig. 1(a). This behavior is quite generic: for 
strong coupling the second order SCPT predicts a nar- 
rower band width than the FC approach, with the error 
of course approaching zero for asymptotically strong cou- 
pling or large phonon frequency, whereas for weaker EP 
coupling, as the convergence of the SCPT deteriorates 
further the band width is overestimated. 

At this point we should discuss the physical content of 
these results. The basic physics that we are observing 
was studied already many years ago by Engelsberg and 
SchrieffcrO for a continuum rather than tight-binding 



model. Considering the single-particle spectral func- 
tion, for weak EP coupling the lowest excitation branch 
at small momentum is a weakly dressed electron. If 
the energy to excite one phonon lies inside the electron 
band, then for arbitrarily weak coupling the electron and 
phonon mix and repel one another near the point where 
they would be degenerate (this may still be quite clearly 
seen in Fig. 1(b), where the dispersion "flattens" at an 
excitation energy near loq). What continues as lowest- 
energy excitation in the single-particle spectral function 
for larger momentum is in fact the phonon, with a small 
admixture of electronic character. The main part of the 
spectral weight disperses upward following the bare elec- 
tron band structure, but is strongly broadened. With in- 
creasing EP interaction the admixture of electronic char- 
acter in the lower energy "flat" part of the dispersion 
at large momentum increases and the peak splits away 
from the continuum, and it begins to be sensible to de- 
note this the polaronic QP band. An interesting new fea- 
ture in the present results is the fact that even in cases 
where the band width of this QP is less than half of the 
bare phonon energy, the QP dispersion may be nearly flat 
at large momentum. In this crossover region low order 
strong-coupling perturbation theory underestimates the 
coherent band width slightly, but more striking is the fact 
that longer range effective hopping terms are not negligi- 
ble, giving rise to a dispersion very different from a n.n. 
cosine band. The SCM approximation on the other hand 
generally overestimates the band width, but qualitatively 
describes the shape of the band quite well. The fact 
that the dispersion is large at small momentum and weak 
at large momentum implies that the mass enhancement 
at the band minimum is generally significantly smaller 
than the exponential prediction of the first order strong- 
coupling perturbation theory m*/m = exp (a 2 ), except 
at very strong coupling. 

We have performed similar calculations also for the 2D 
square lattice. In this case the application of the FC ap- 
proach requires separate determination of h^ n > for n-site 
clusters of all inequivalent topologies, and in the applica- 
tion of h all possible shapes that topologically equivalent 
embedded subclusters may have must be allowed for. At 
the present level of approximation this is still quite sim- 
ple, however, as there are still only three inequivalent 
topologies for connected clusters of 5 sites. The limiting 
factor is the size of system for which sufficiently accurate 
results can be obtained for the eigenvalues and eigenvec- 
tors of the Holstein model, and not yet the number of 
graphs which need be considered. To a good approxima- 
tion the results for 2D are similar to those for ID, af- 
ter renormalization of parameters to take account of the 
doubling of the bare electronic band width in 2D when 
compared to ID. The convergence of the FC approach 
in 2D for u>a/t = 2 and a given coupling is similar to 
that for ujo/t — 1 and the same coupling in ID. In Fig. 
2 we present the QP band structure in 2D for EP cou- 
pling g/t = 2.75 and phonon frequency u>o/t — 1.5 (a) 
and ojo/t = 2.0 (b). In both cases, from a comparison of 
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the 4- and 5-site FC approximations it is again clear that 
the polaron size is still slightly larger than the range of 
effective interactions presently included, but that further 
corrections will be quite small. For the larger phonon fre- 
quency of Fig. 2(b) the discrepancy between second order 
SCPT and the FC is quite small, consistent with good 
convergence of the SCPT for large ojo- Since the effective 
coupling a is not very large in this case, also the SCM 
result is reasonable. For the smaller phonon frequency 
of Fig. 2(a) however, a is significantly larger, and the 
band width predicted by the FC approximation is much 
smaller than the weak-coupling SCM result. Once again 
the first and second order SCPT err in the opposite di- 
rection compared to the SCM result, and underestimate 
the band width by more than a factor of two. 




(0,0) 



Fig. 2. Polaron dispersion for 2D system with g/t = 2.75 
and: (a) u>o/t = 1.5; (b) ujo/t = 2.0. Solid lines are FC ap- 
proximation to h^ 5 \ dotted lines FC to h^K Dashed lines are 
first order strong-coupling perturbation theory, long dashed 
lines are second order strong-coupling perturbation theory, 
and dot-dashed lines weak-coupling self-consistent Migdal ap- 
proximation. 

However, already in the second order SCPT correction 
the tendency for the shape of the dispersion to "flatten" 
at large momentum when compared with the n.n. co- 
sine (first order SCPT) is apparent. One should again 
stress that the extended flat region at large momentum 
for such parameters is not trivial: the FC band width is a 
factor of three smaller than the weak-coupling SCM pre- 
diction w u>o, where the flat region is clearly determined 
by the "collision" of the bare electronic dispersion with 
the phonon. 

We have presented an approach to the polaron problem 
in the Holstein model which is equivalent to a resumma- 
tion of some terms in strong-coupling perturbation the- 
ory to all orders. For phonon frequencies uJo/t greater 
than 1 (2) in ID (2D) this approach is sufficiently con- 



vergent that good agreement can be achieved with the 
SCM weak-coupling perturbation theory in the crossover 
region. The flat band at large momentum, expected for 
weak EP coupling due to the hybridization of the elec- 
tron and phonon degrees of freedom when the phonon 
frequency lies inside the electron band, has been shown to 
survive to surprisingly strong EP coupling. For too/t ~ 1 
and intermediate coupling the low-order strong-coupling 
perturbation theory has been shown to underestimate the 
QP band width, and to overestimate the mass enhance- 
ment by an even larger factor due to the importance 
of longer range hopping terms in the effective polaron 
model. 
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